Donor–Acceptor–Donor 1H-Benzo[d]imidazole Derivatives as Optical Waveguides

A new series of donor–acceptor–donor (D–A–D) structures derived from arylethynyl 1H-benzo[d]imidazole was synthesized and processed into single crystals with the goal of testing such crystals’ ability to act as optical waveguides. Some crystals displayed luminescence in the 550–600 nm range and optical waveguiding behavior with optical loss coefficients around 10−2 dB/μm, which indicated a notable light transport. The crystalline structure, confirmed by X-ray diffraction, contains internal channels that are important for light propagation, as we previously reported. The combination of a 1D assembly, a single crystal structure, and notable light emission properties with low losses from self-absorption made 1H-benzo[d]imidazole derivatives appealing compounds for optical waveguide applications.


Introduction
In recent years, photonic integrated circuits (PICs) have become of paramount importance in many applications, such as optical signal processing, biomedical applications, and optical computing. Although the most common platforms for PICs are now based on inorganic bulk semiconductors such as InP or Si, alternatives such as 2D semiconductors [1], colloidal quantum dots [2], hybrid perovskites [3] and organic crystals [4] are the subjects of extensive research.
Optical waveguide is a pivotal and fundamental component of PICs. This physical structure can conduct and propagate light efficiently with minimal loss. The two factors that influence the behavior of the optical waveguide are the total internal reflection and the fact that the medium must have an index of refraction higher than the external medium. Thus, when the incident light angle is higher than the critical angle, the light is confined and propagated along the structure, following the Snell law [5,6].
Organic crystals possess unquestionable advantages as light emitting and transporting media, including their capability to guide light in the visible region giving rise to large operational bandwidths, low transmission losses, simple and inexpensive processing into micro-/nanofibers of interest for miniaturization, high refractive indices, well-defined crystal faces with few grain boundaries and surface defects, and better heat dissipation [7][8][9]. Owing to their large chemical versatility, the waveguide geometry and transmission properties can be tailored through chemical functionalization of the molecular building blocks. For these reasons, in recent years there has been a large increase in the use of organic crystals as optical waveguides [10][11][12][13][14][15][16].
Integration of organic crystals in PICs remains a challenge. Nevertheless, extensive research is being carried out in this field to achieve such integration. For example, a PIC was achieved by using two different organic crystals in which, upon receiving the signal input, the circuit executed direction-specific optical outputs [17]. Furthermore, it is now possible to identify various real applications in which optical waveguides have already been successfully implemented, including highlighting the novel 5G technology with the employment of organic flexible crystals [18], use as media in optical communications [19], application in the construction of optical waveguide multimodal sensors for measuring human temperature [20], or as materials that possess the ability to convert light into mechanical work [21], demonstrating the relevance of this type of physical structure in the present and in the near future.
We recently described some examples of organic molecules with D-A-D architectures that show optical waveguiding behavior. They consist of an electron donor (D) and an electron acceptor (A) connected by a π-conjugated bridge [22], which facilitates intramolecular charge transfer (ICT) [23] and allows modulation of the energy bandgap, giving rise to luminescent optical waveguides that emit over a wide spectral range [22].
In this report, we extend this concept to the 1H-benzo [d]imidazole moiety, an interesting core, as it displays a strong acceptor character due to the presence of a C=N bond, which can be easily modulated with the introduction of different groups in positions 4 and 7. In addition, the NH group facilitates crystalline self-assembling via intermolecular hydrogen bonding, which is required for obtaining fibers or crystals with optical waveguiding behavior.
Despite these interesting properties, 1H-benzo [d]imidazole is still an unexplored moiety in material science, and only some examples of 1H-benzo[d]imidazole derivatives can be found, with applications as electrochromic materials [24], in organic electronics as photoactive components in organic solar cells [25] or OLEDs [26], and in gels [27].
Two electron-attracting trifluoromethyl groups were introduced in the 2-aryl group in order to increase the acceptor character of the benzimidazole moiety. In the same way, three different donor groups were introduced in positions 4 and 7 using alkynyl groups as linkers, as alkynyl groups facilitate the formation of supramolecular structures due to their ability to form CH-π interactions [22,28]. It should be noted that, to date, no 1Hbenzo[d]imidazole-derived crystals with optical waveguide properties have been reported.

Preparation of 1H-Benzo[d]imidazole Derivatives 1
1H-benzo[d]imidazole derivatives 1 were synthesized in two steps, as shown in Scheme 1. Thus, a mixture of 3,6-dibromobenzene-1,2-diamine (2) and 3,5- Two electron-attracting trifluoromethyl groups were introduced in the 2-aryl group in order to increase the acceptor character of the benzimidazole moiety. In the same way, three different donor groups were introduced in positions 4 and 7 using alkynyl groups as linkers, as alkynyl groups facilitate the formation of supramolecular structures due to their ability to form CH-π interactions [22,28]. It should be noted that, to date, no 1Hbenzo[d]imidazole-derived crystals with optical waveguide properties have been reported.

Preparation of 1H-Benzo[d]imidazole Derivatives 1
1H-benzo[d]imidazole derivatives 1 were synthesized in two steps, as shown in Scheme 1. Thus, a mixture of 3,6-dibromobenzene-1,2-diamine (2) and 3,5-bis(trifluoromethyl) benzaldehyde (3) in the presence of p-toluenesulfonic acid resulted in the dibromo derivative 4 in 84% yield. Synthesis of the final bisarylethynyl-1H-benzo[d]imidazoles 1a-c was attempted by a microwave-assisted Sonogashira C-C cross-coupling reaction between 4 and the corresponding alkynyl derivative 5, as we previously described [22]. However, unreacted dibromo benzimidazole 4 was recovered in a complex mixture. Recently, it was stated that the presence of pyridinic or pyrrolic nitrogens and their different Cu(I) coordination may influence the activation of the triple bond [29]. Then, the basicity of the pyridinic nitrogen of the benzimidazole could also justify the lack of reactivity in the Sonogashira reaction. Derivative 1 was obtained in good yields (53-87%) using an alternative microwave-assisted Stille reaction between 4 and the corresponding arylethynyllstannane 6.
Molecules 2023, 28, x FOR PEER REVIEW 3 of 15 bis(trifluoromethyl)benzaldehyde (3) in the presence of p-toluenesulfonic acid resulted in the dibromo derivative 4 in 84% yield. Synthesis of the final bisarylethynyl-1Hbenzo[d]imidazoles 1a-c was attempted by a microwave-assisted Sonogashira C-C crosscoupling reaction between 4 and the corresponding alkynyl derivative 5, as we previously described [22]. However, unreacted dibromo benzimidazole 4 was recovered in a complex mixture. Recently, it was stated that the presence of pyridinic or pyrrolic nitrogens and their different Cu(I) coordination may influence the activation of the triple bond [29]. Then, the basicity of the pyridinic nitrogen of the benzimidazole could also justify the lack of reactivity in the Sonogashira reaction. Derivative 1 was obtained in good yields (53-87%) using an alternative microwave-assisted Stille reaction between 4 and the corresponding arylethynyllstannane 6.

Theoretical Calculations
A deeper understanding of the non-reactivity shown by 4 in Sonogashira reactions, as compared to analogous benzotriazole systems, has been achieved by analysis of the aromaticity of the rings within the DFT framework at B3LYP-D3/6-31G(d) level via the nucleus-independent chemical shift (NICS) method developed by Schleyer et al. [30]. NICS has been widely used as an aromaticity probe by computing the absolute magnetic shielding at the point of interest [31]. In order to have a complete overview, we calculated the NICS values at the ring-critical points and at 1 Å above as NICS(0) and NICS(1), respectively. The results obtained for the three aromatic rings of 4 are shown in Figure 2.
The NICS(0) and NICS(1) values obtained for ring C (−11.21 and −10.93 ppm, respectively) were more negative than the values described for analogous dibrominated benzotriazole (NICS(0) and NICS(1) values of -8.92 and -9.47 ppm, respectively) previously employed by our group [22]. Therefore, the NICS probe indicated a higher aromaticity of the dibrominated core in 4 compared to analogous dibrominated benzotriazole, in line with the lower reactivity of the former toward the Sonogashira coupling reaction (i.e., harder conditions are required for 1H-benzo[d]imidazole core, compared to benzotriazole analogous).

Theoretical Calculations
A deeper understanding of the non-reactivity shown by 4 in Sonogashira reactions, as compared to analogous benzotriazole systems, has been achieved by analysis of the aromaticity of the rings within the DFT framework at B3LYP-D3/6-31G(d) level via the nucleus-independent chemical shift (NICS) method developed by Schleyer et al. [30]. NICS has been widely used as an aromaticity probe by computing the absolute magnetic shielding at the point of interest [31]. In order to have a complete overview, we calculated the NICS values at the ring-critical points and at 1 Å above as NICS(0) and NICS(1), respectively. The results obtained for the three aromatic rings of 4 are shown in Figure 2.
The NICS(0) and NICS(1) values obtained for ring C (−11.21 and −10.93 ppm, respectively) were more negative than the values described for analogous dibrominated benzotriazole (NICS(0) and NICS(1) values of -8.92 and -9.47 ppm, respectively) previously employed by our group [22]. Therefore, the NICS probe indicated a higher aromaticity of the dibrominated core in 4 compared to analogous dibrominated benzotriazole, in line with the lower reactivity of the former toward the Sonogashira coupling reaction (i.e., harder conditions are required for 1H-benzo[d]imidazole core, compared to benzotriazole analogous).
On the other hand, with the aim of analyzing the molecular properties of the synthesized compounds, the minimum-energy optimized structures at B3LYP/6-31G(d,p) theoretical level were calculated. In all cases, the highest occupied molecular orbital (HOMO) was located mainly in the horizontal branch of the molecules, while the lowest unoccupied molecular orbital (LUMO) was located mainly in the vertical branch, with an overlap between both orbitals on the central acceptor core, which facilitates intramolecular charge transfer (ICT) from the electron-donating to the electron-withdrawing units. The calculated HOMO-LUMO energies, HOMO-LUMO gaps and topologies are shown in Figure 3. On the other hand, with the aim of analyzing the molecular properties of the synthesized compounds, the minimum-energy optimized structures at B3LYP/6-31G(d,p) theoretical level were calculated. In all cases, the highest occupied molecular orbital (HOMO) was located mainly in the horizontal branch of the molecules, while the lowest unoccupied molecular orbital (LUMO) was located mainly in the vertical branch, with an overlap between both orbitals on the central acceptor core, which facilitates intramolecular charge transfer (ICT) from the electron-donating to the electron-withdrawing units. The calculated HOMO-LUMO energies, HOMO-LUMO gaps and topologies are shown in Figure 3.   On the other hand, with the aim of analyzing the molecular properties of the synthesized compounds, the minimum-energy optimized structures at B3LYP/6-31G(d,p) theoretical level were calculated. In all cases, the highest occupied molecular orbital (HOMO) was located mainly in the horizontal branch of the molecules, while the lowest unoccupied molecular orbital (LUMO) was located mainly in the vertical branch, with an overlap between both orbitals on the central acceptor core, which facilitates intramolecular charge transfer (ICT) from the electron-donating to the electron-withdrawing units. The calculated HOMO-LUMO energies, HOMO-LUMO gaps and topologies are shown in Figure 3.  Modification of the peripheral donor groups produced variations in the energy values of the frontier molecular orbitals, so that the HOMO-LUMO gap could be modulated by introducing different electron-donating groups. Thus, compound 1a, which possesses weaker donor groups, showed a higher HOMO-LUMO gap value than derivatives 1b and 1c, with stronger donor groups and more effective π-conjugation.

Photophysical Data
The UV-vis and fluorescence spectra of derivatives 1, recorded in 10 −5 M CHCl 3 solutions, are shown in Figure 4 and the photophysical values are summarized in Table 1.
The UV-vis and fluorescence spectra of derivatives 1, recorded in 10 M CHCl3 solutions, are shown in Figure 4 and the photophysical values are summarized in Table 1. 1a-c shows an absorption band ascribed to a π-π* transition (284-316 nm) and a dominant band at low energy (342-384 nm), ascribed to an intramolecular charge transfer (ICT), as supported by computational studies ( Figure S1 and Table S1). Note that the theoretical results were in good agreement with the experimental results. Displacement of the absorption bands (λmax) followed the 1a < 1b < 1c trend, which was related to the increase in the donor character of the peripheral groups. Thus, phenyl derivative (1a) showed an absorption band shifted to the blue, as a consequence of the weaker electrondonor character of the phenyl group, while 1b (trimethoxyphenyl) and 1c (benzothiophenyl), with higher donor character, resulted in a bathochromic shift of the maxima.  The fluorescence bands of 1a-c were between 429-459 nm and followed the same trend as absorption. In contrast to the absorption maxima, 1b showed a slightly larger shift of the emission band in comparison with 1c.
In general, all compounds displayed moderate-to high-fluorescence quantum yields ranging from 0.55 to 0.78.   1a-c shows an absorption band ascribed to a π-π* transition (284-316 nm) and a dominant band at low energy (342-384 nm), ascribed to an intramolecular charge transfer (ICT), as supported by computational studies ( Figure S1 and Table S1). Note that the theoretical results were in good agreement with the experimental results. Displacement of the absorption bands (λ max ) followed the 1a < 1b < 1c trend, which was related to the increase in the donor character of the peripheral groups. Thus, phenyl derivative (1a) showed an absorption band shifted to the blue, as a consequence of the weaker electron-donor character of the phenyl group, while 1b (trimethoxyphenyl) and 1c (benzothiophenyl), with higher donor character, resulted in a bathochromic shift of the maxima.
The fluorescence bands of 1a-c were between 429-459 nm and followed the same trend as absorption. In contrast to the absorption maxima, 1b showed a slightly larger shift of the emission band in comparison with 1c.
In general, all compounds displayed moderate-to high-fluorescence quantum yields ranging from 0.55 to 0.78.

Molecular Self-Assembly
The generation of supramolecular aggregates was carried out by the slow diffusion technique. Thus, a vial containing a dilute solution (1 mg) of 1a-c in 1 mL of tetrahydrofuran or chloroform, as good solvents, was gently introduced into a second vessel containing a poor solvent, such us hexane, acetonitrile, methanol, or ethanol. This second vessel was closed and, after seven days, the formation of supramolecular structures was observed.
1a did not lead to the formation of supramolecular structures in any of the mixtures used. In contrast, 1b resulted in the formation of well-defined needle-shaped crystals furan or chloroform, as good solvents, was gently introduced into a second vessel contain-ing a poor solvent, such us hexane, acetonitrile, methanol, or ethanol. This second vessel was closed and, after seven days, the formation of supramolecular structures was observed.
1a did not lead to the formation of supramolecular structures in any of the mixtures used. In contrast, 1b resulted in the formation of well-defined needle-shaped crystals using different solvent mixtures (Figure 5a-c), whereas 1c only crystallized in THF/ethanol (Figure 5d).
The thicknesses of the structures ranged from 5 to 20 µm and lengths ranged from tens to hundreds µm, in some cases reaching more than 1 mm in length.

Optical Waveguide Properties
To elucidate their optical waveguide properties, crystals of 1b and 1c were studied by fluorescence microscopy in the experimental setup shown in Figure 6. Thus, the crystals were excited in the fiber body with a photoexcitation centered at 365 nm. The thicknesses of the structures ranged from 5 to 20 µm and lengths ranged from tens to hundreds µm, in some cases reaching more than 1 mm in length.

Optical Waveguide Properties
To elucidate their optical waveguide properties, crystals of 1b and 1c were studied by fluorescence microscopy in the experimental setup shown in Figure 6. Thus, the crystals were excited in the fiber body with a photoexcitation centered at 365 nm. Crystals of 1b obtained in CHCl3/CH3CN (Figure 7a) and 1c in THF/ethanol ( Figure  7b) showed bright yellow and green fluorescence, respectively. The PL spectra of 1b (Figure 7c) and 1c (Figure 7d) crystals were rather broad, with maxima at 586 nm and 564 nm, respectively. The bright emission observed at the fiber tip was motivated by outcoupled photons, due to the refractive index mismatch at the crystal/air interface, and confirmed that luminescence was efficiently waveguided along the fiber axis.  Crystals of 1b obtained in CHCl3/CH3CN (Figure 7a) and 1c in THF/ethanol ( Figure  7b) showed bright yellow and green fluorescence, respectively. The PL spectra of 1b (Figure 7c) and 1c (Figure 7d) crystals were rather broad, with maxima at 586 nm and 564 nm, respectively. The bright emission observed at the fiber tip was motivated by outcoupled photons, due to the refractive index mismatch at the crystal/air interface, and confirmed that luminescence was efficiently waveguided along the fiber axis. In both compounds, a bathochromic shift of the emission bands was observed in the solid state with respect to solution motivated by the enhanced dielectric constant in the solid medium, as well as by intermolecular interactions [32]. Moreover, in both solid state and in solution, 1b showed a larger red shift than 1c, indicating the influence of the donor group in both states.
Next, we evaluated the optical loss coefficient (OLC) of 1c fibers only, as 1b fibers were too small to be measured with our experimental setup capabilities.
To estimate the OLC of 1c, a representative crystal of length ≈ 450 µm was studied. For this purpose, the crystal was excited with a 355 nm laser at different positions along the longitudinal axis and the resulting emission was detected at the fiber tip. As the light propagation distance increased, the spectra showed a progressive decrease in photoluminescence (PL) intensity caused by self-absorption, following the Beer-Lambert law, enabling the obtaining of the absorption coefficient (α) in µm −1 .
bling the obtaining of the absorption coefficient (α) in µm −1 .

X-ray Studies: Structure-Property Relationships
We have reported some structural requirements for supramolecular aggregates to behave as optical waveguides, adding a new dimension to the study of organic optical waveguides [34]. In this sense, the presence of internal channels should favor the propagation of light through the crystal. Furthermore, recent literature showed that the presence of microcavities within the crystals promotes the strong coupling between excitons and polaritons, resulting in the formation of exciton polaritons (Eps) that are responsible for light propagation [36].
Consequently, in order to further deepen the knowledge of the structure/property relationship in optical waveguides, an X-ray study of the 1b,c crystals structures was carried out.
Unfortunately, only the crystal structure of 1b was successfully resolved by single crystal X-ray diffraction. Compound 1b crystallized in the Pī space group and two molecules of 1b and two molecules of acetonitrile used for crystallization completed the packing of the unit cell.
The conformation of 1b was non-planar, as the peripheral rings were rotated (7.3° for the A-ring, 70.3° for the B-ring and 15.1° for the C-ring) with respect to the 1H-

X-ray Studies: Structure-Property Relationships
We have reported some structural requirements for supramolecular aggregates to behave as optical waveguides, adding a new dimension to the study of organic optical waveguides [34]. In this sense, the presence of internal channels should favor the propagation of light through the crystal. Furthermore, recent literature showed that the presence of microcavities within the crystals promotes the strong coupling between excitons and polaritons, resulting in the formation of exciton polaritons (Eps) that are responsible for light propagation [36].
Consequently, in order to further deepen the knowledge of the structure/property relationship in optical waveguides, an X-ray study of the 1b,c crystals structures was carried out.
Unfortunately, only the crystal structure of 1b was successfully resolved by single crystal X-ray diffraction. Compound 1b crystallized in the Pī space group and two molecules of 1b and two molecules of acetonitrile used for crystallization completed the packing of the unit cell.
The conformation of 1b was non-planar, as the peripheral rings were rotated (7.3 • for the A-ring, 70.3 • for the B-ring and 15.1 • for the C-ring) with respect to the 1Hbenzo[d]imidazole core (Figure 9a). Therefore, 1b molecules were arranged coplanar and antiparallel at 3.5 Å from each other (Figure 9b). This arrangement was supported by hydrogen bonds between the pyridinic nitrogen and H ortho (2.49 Å) and with the methoxy groups of the B-ring (2.72 Å). Moreover, the methoxy groups of the B-ring interacted via hydrogen bond with the homologous groups of the A-ring (2.40 Å), as well as with the -CF 3 groups (2.60 Å) (Figure 9c).
The formation of molecular layers was based on hydrogen bonds, where acetonitrile acted as a linker between 1b molecules by interacting with the -NH group (2.04 Å) and the H ortho (2.46 Å) and methoxy groups of the A-ring (2.60 Å). C-H···π interactions between the triple bonds and the C-ring (3.39 Å) allowed the stacking of the molecular layers (Figure 9d).
In the three-dimensional arrangement of 1b, two types of channels (A and B) extending in the direction of the a axis were observed. As shown in Figure 10, channels A were located inside the unit cell and channel B was located on the edges of the unit cell. oxy groups of the B-ring interacted via hydrogen bond with the homologous groups of the A-ring (2.40 Å), as well as with the -CF3 groups (2.60 Å) (Figure 9c).
The formation of molecular layers was based on hydrogen bonds, where acetonitrile acted as a linker between 1b molecules by interacting with the -NH group (2.04 Å) and the Hortho (2.46 Å) and methoxy groups of the A-ring (2.60 Å). C-H···π interactions between the triple bonds and the C-ring (3.39 Å) allowed the stacking of the molecular layers (Figure 9d). In the three-dimensional arrangement of 1b, two types of channels (A and B) extending in the direction of the a axis were observed. As shown in Figure 10, channels A were located inside the unit cell and channel B was located on the edges of the unit cell. The formation of channel A can be attributed to hydrogen-bonding interactions between the two molecules that complete the packing of the unit cell ( Figure S2, red lines). On the other hand, formation of channel B was exclusively due to the presence of C-H … π interactions between a C26 methoxy group and 3,5-bis(trifluoromethyl)phenyl ring, which occurred along the b axis ( Figure S2, blue lines).
The extension of the channels along the a axis and perpendicularly to the longest axis of growth (c axis) occurred due to the π … π interactions between 3,5-bis(trifluoromethyl)phenyl and 3,4,5-(trimethoxy)phenyl rings, with a distance of 4.05 Å between their planes ( Figure S3).
These results confirmed that the presence of internal channels located in the direction perpendicular to the longest axis of growth allows the propagation of the light, supporting our previously reported hypothesis [34]. The formation of channel A can be attributed to hydrogen-bonding interactions between the two molecules that complete the packing of the unit cell ( Figure S2, red lines). On the other hand, formation of channel B was exclusively due to the presence of C-H . . . π interactions between a C26 methoxy group and 3,5-bis(trifluoromethyl)phenyl ring, which occurred along the b axis ( Figure S2, blue lines).
The extension of the channels along the a axis and perpendicularly to the longest axis of growth (c axis) occurred due to the π . . . π interactions between 3,5-bis(trifluoromethyl)phenyl and 3,4,5-(trimethoxy)phenyl rings, with a distance of 4.05 Å between their planes ( Figure S3).
These results confirmed that the presence of internal channels located in the direction perpendicular to the longest axis of growth allows the propagation of the light, supporting our previously reported hypothesis [34].

Materials and Reagents
The reagents needed for the synthesis of all the products described in this article that were commercially available were used with any kind of prior purification. All the air-sensitive reactions, such as the formation of the stannanes or the Still cross coupling reactions, were performed employing an argon atmosphere, making three vacuum-argon cycles. Synthesis in which microwave irradiation was the energy source was performed in a microwave Discover ® (CEM, Matthews, NC, USA).
NMR spectra ( 1 H and 13 C) were performed on Bruker Advance Neo NMR spectrometers working at 399.77 and 500.16 MHz for 1 H and 100.53 and 125.75 MHz for 13 C, respectively. The spectra were taken at 298 K in all the cases, employing partially deuterated solvents as the internal reference. The coupling constants (J) were described in Hertz (Hz) and chemical shifts (δ) in ppm. Signals' multiplicities are noted as s = singlet, d = doublet, t = triplet, q = quadruplet, quint = quintuplet, sext = sextet, sept = septet, and m = multiplet.
For the photophysical studies, UV-visible studies were performed in a Jasco V-750 spectrophotometer (JASCO-Spain, Madrid, Spain) and the fluorescence studies were performed in a Jasco FP-8300 spectrofluorometer (JASCO-Spain, Madrid, Spain). The absorption and emission spectra shown in the main text were obtained from 10 −5 M chloroform (with high spectroscopic grade) solutions at room temperature employing standard quartz cells with a width of 1 cm.
In order to determine the melting points, a Büchi model M-569 melting point meter was employed.
Regarding mass spectroscopy, MALDI-TOF mass spectra were employed using a Bruker Autoflex II TOF/TOF spectrometer with the employment of a matrix of dithranol, employing the positive mode. All the samples, which were co-crystallized with the matrix on the probe, were then ionized via a nitrogen laser pulse of 337 nm and accelerated under 20 Kv with time-delayed extraction before entering the time-of-flight mass spectrometer. Matrix (10 mg/mL) and sample (1 mg/mL) were separately dissolved in tetrahydrofuran (they are soluble in this solvent) and mixed in a different matrix/sample ratio range from 100:1 to 50:1. As is typical in these kinds of measurements, a 5 mL mixture of matrix and sample was applied to a MALDI-TOF MS probe and dried later under air. In addition, external calibration was carried out with the employment of Peptide Calibration Standard II (covered mass range: 700-3200 Da) from Care (Bruker, Billerica, MA, USA). The applied peak (m/z determination) detection method was the threshold centroid at 50% height of the maximum peak.
SEM images to determine the morphology of the supramolecular aggregates were obtained from a SEM Zeiss Gemini SEM 500 (Carl Zeiss AG, Oberkochen, Germany), which operated at 3 kV. The different crystals of all the derivatives, obtained by the slow diffusion technique, were deposited onto a glass substrate and the remaining solvent was slowly evaporated, without any further treatment.
X-ray diffraction data were extracted from a Bruker X8 APEX II CCD-based diffractometer (Bruker, Billerica, MA, USA). Data were processed using SAINT and, after that, an absorption correction was performed with the program SADABS. The structures were solved mixing direct methods and difference Fourier syntheses and, finally, refined by fullmatrix least-squares on F2 employing the WINGX version 2021.3 and OLEX2-1.5 software packages. The rest of the atoms, except hydrogen, were refined with anisotropic thermal parameters. Hydrogen atoms were placed using a "riding model" and were included in the refinement at calculated positions.
To study the optical waveguiding behavior, PL microscopy images were taken from a Nikon Eclipse Ti inverted microscope (Nikkon Corporation, Tokyo, Japan) with dry objectives (100× N.A. 0.8 and 20× N.A. 0.45) coupled to a Shamrock spectrometer from Andor Technology (Andor, Belfast, UK) with a thermoelectrically cooled Newton EM (Andor) CCD. The excitation was obtained taking into account an appropriate filtering of the lines from a Xe lamp.

Procedure for the Synthesis of Arylakynylstannanes 6
Arylethynylstannane 6a is commercially available and was used without further purification. From 5b-c, the arylethynylstannanes 6b-c were synthesized, as previously reported [38].
General procedure: To a solution of the corresponding acetylene derivative (5b-c) (2 mmol) in dry THF (−78 • C) and under inert atmosphere, n-BuLi (1.87 mL, 3 mmol) was added slowly dropwise. Then, the mixture was stirred for 30 min. and, at that moment, Bu 3 SnCl (0.64 mL, 2.4 mmol) was added dropwise. The reaction mixture was allowed to warm until reaching room temperature and stirred again for 2 h. Then, to the crude reaction, 30 mL of water was added, neutralized with HCl 1M, and the organic layers were extracted with dichloromethane (3 × 30 mL). The organic layers were dried with MgSO 4 , the solvent was removed in vacuo, and the products were employed in the next step without further purification.  A mixture of 4 (0.12 g, 0.25 mmol), the corresponding ethynylstannane 6 (0.62 mmol), PdCl 2 (PPh 3 ) 2 (0.007 g, 0.01 mmol) and LiCl (0.03 g, 0.75 mmol) was charged under inert atmosphere to a microwave vessel. CH 3 CN (1 mL) was added, the vessel was closed and irradiated at 110 • C for 20 min. The crude reaction was purified by column chromatography on silica gel, eluting with hexane/ethyl acetate to afford the pure products 1a-c.

Conclusions
In this article, we described a new series of different symmetrical D-A-D 4,7bis(arylethynyl)-1H-benzo[d]imidazole derivatives, introducing different donors and acceptors into the structure to modulate the energy value of the HOMO-LUMO gap.
SEM images of the assembled crystals indicated that 1b and 1c derivatives formed needle-shaped crystals, suitable for application as optical waveguides.
The PL microscopy studies revealed that 1b and 1c showed strong yellow-to-green luminescence with PL maxima around 586 and 564 nm, respectively. In addition, 1c showed a low OLC, exhibiting an efficient light propagation. Thus, compounds 1b and 1c are the first reported 1H-benzo[d]imidazole derivatives to be applicable as optical waveguides.
Finally, the excellent morphology of an aggregate of 1b allowed its study by X-ray diffraction. An in-depth analysis of the crystal structure verified that the existence of internal channels inside the crystal allows the transmission of light for the optical waveguiding behavior. These results are in agreement with previous X-ray diffraction studies on structures reported by our research group, confirming the first approximation of the structure/optical waveguide properties relationship.